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Abstract 

In the paper jl3) , Stanley gives formulas for the number of plane partitions in each of ten symmetry 
classes. This paper together with the two papers |lj and JliJ completes the project of proving all ten 
formulas. 

We enumerate cyclically symmetric, self-complementary plane partitions. We first convert plane 
partitions to tilings of a hexagon in the plane by rhombuses, or equivalently to matchings in a certain 
planar graph. We can then use the permanent-determinant method or a variant, the Hafhian-Pfafnan 
method, to obtain the answer as the determinant or PfafHan of a matrix in each of the ten cases. We 
row-reduce the resulting matrix in the case under consideration to prove the formula. A similar row- 
reduction process can be carried out in many of the other cases, and we analyze three other symmetry 
classes of plane partitions for comparison. 

1 Introduction 

In the paper |p"3f , Richard Stanley describes ten symmetry classes of plane partitions and gives formulas for 
their enumeration. We give a brief summary of the conjectures and results mentioned in the paper: 

A plane partition is a finite union 7r of unit cubes in the positive octant of M 3 such that for each cube C 
in 7T, there is either another cube or a wall of the octant below C, and the same is true behind and to the 
left of C . (Sometimes this arrangement is called the Ferrer's diagram of a plane partition.) Alternatively, a 
plane partition is a finite subset of Z + x Z + x Z + which is an order ideal with respect to the partial ordering 
(x + 1, y, z), (x, y + 1, z), (x, y, 2 + 1) > (x, y, z) for (x, y, z) € Z^_. An example of an ordinary but useful plane 
partition is a box B(a, b, c) = [0, a] x [0, b] x [0, c]. We define three actions on the set of pairs (tt, B(a, b, c)), 
where 7r C B(a, 6, c): 

1. We define t(tt), the transpose of tt, to be {(x, y, z)| (y, x, z) G 7r}, and t(tt, B(a, b, c)) to be (r(7r), B(b, a, c)). 

2. We define p(tt), the rotation of tt, to be {(x, y, z) | (y, z, x) £ tt}, and p(n, B(a, b, c)) to be (p(ir), B(c, a, b)). 

3. We define k(tt, B(a,b,c)) = (ir 1 , B(a, b, c)), where ir' is the complement of tt in B(a,b,c) with all three 

coordinates reversed. Whenever a box B(a,b,c) has been chosen, we may write for tt' by abuse 
of notation. 

Let T be the group generated by r, p, and k. For each subgroup G of T, we define Nc(a,b,c) to 
be the number of pairs (tt, B(a, b, c)) fixed by G (in general, Nc{a,b,c) is zero unless B(a,b,c) is fixed 
by G). There are ten inequivalent choices Gi, . . . , Gio for G; we set Ni(a,b,c) — Na i (a,b,c). Table [I] 
gives the definition of each of the ten groups, together with the corresponding formula for iVj. In the 
table, H(n) is the hyperfactorial function and Hk(n) is a staggered hyperfactorial, defined by Hk{n) — 
(n — k)l(n — 2k)\(n — 3fc)! . . ., H(n) = Hi(n), and Fk(n) — n(n — k)(n — 2k) ... is a staggered factorial. The 

symbol = indicates that the equality was conjectural before this paper. 

In this paper, we will obtain an expression for each of the ten numbers N%, . . . , A^o as the determinant 
or the Pfaffian of a matrix. We will then row- and column-reduce the matrices to prove the formulas in 
Table | in cases 1, 3, (part of) 5, and 9. The author has also found a way to evaluate the determinants 
of the matrices for cases 2, 6, 7, and 8, but remains confounded by cases 4 and 10. The author feels that 
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Group 


Formula 


Gi = (e) 
G-2 = (r) 
G 3 = (p) 
G 4 = (t, p) 


Ni{a,b lC ) 
N 2 (a,a,b) 
N 3 (a, a, a) 
Ni{a, a, a) 


H(a+fc+c)_ff(a)_H"(6)//(c) 
ff(a+b)_ff(a+c)if(b+c) 

H 2 (2a+b+l)_ff(a)_f/ 2 (b) 
Ht (2a+l)H(a+b) 

H 3 (3o+2)H(a) 
H(2a)F 3 (3a-2) 
_ // 2 (a)_ff 6 (3a+5) 
H 2 (2a+l)_F 6 (3a-2) 


G 5 = (k) 


N 5 (2a,2b,2c) 


iVi(a,6,c) 2 




N 5 (2a,2b, 2c + 1) 


Ni(a,b,c)Ni(a,b,c+ 1) 




N 5 (2a+ 1,2b + 1,2c) 


= iVi(a + 1,6, c)7Vi(a,6+ l,c) 


G 6 = (*r> 


iV 6 (a,a,26) 


H 2 (2fc+l)H 2 (26+2a)H(a) 
H(2b+a)ff 2 (2a) 


G 7 = <k,t> 


N 7 (2a,2a,26) 


= iVi(a, a, 6) 




7V 7 (2a + 1,2a + 1,26) 


7Vi(a,a+l,6) 


G 8 = (kt, p) 

Gg = (K, p) 


N 8 (2a,2a,2a) 
W 9 (2a,2a,2a) 


F 3 (3a-2)H 6 (6a)H 2 (2a) 
~ H 4 (4a+l)ff 4 (4a) 
! H 3 (3a+l) 2 H(a) 2 
ff(2a) 2 


Gio = (k,t,p) 


iYi (2a,2a,2a) 


ff 3 (3a+l)H(a) 
ff(2a) 



Table 1: Enumeration formulas for symmetric plane partitions 



the permanent-determinant method is a good first step towards a unified treatment of all ten cases, but the 
entirety of such a treatment remains elusive. 

The proofs of cases 4 and 10 are more recent than the results of this paper and will appear in |l4| and . 
These three results conclude the problem of enumeration of plane partitions in different symmetry classes, 
although the problem of g-enumeration (see section remains open in case 4. Table || gives an account of 
who first conjectured and who first proved the formula in each symmetry class. In particular, the formula in 
case 9, whose proof is the main news of this paper, i Table |3| gives the same information for g-enumeration. 
The tables are largely compiled from jl3| . They list the number of the symmetry class as given by Stanley, 
the name in the convention of Stembridge, and the generators in the notation of this paper. 

Our plan is to combine two old ideas about combinatorial tilings. The first idea has probably been known 
since antiquity, and is illustrated in || and Jl^] and discussed in great generality in fL6f . It is that a plane 
partition in a box is equivalent to a tiling by certain rhombuses of a certain hexagon. To convert a plane 
partition to a tiling, we view the partition as a union of unit cubes in an a x 6 x c box, and we simply "draw 



Stanley 


Stembridge 


Group 


First formulator 


First enumerator 


1 


P 


<e> 


MacMahon 


MacMahon 


2 


S 


<r) 


MacMahon 


Andrews, Gordon, Macdonald 


3 


CS 


(P) 


Macdonald 


Andrews 


4 


TS 


(P,r) 


Stanley 


Stembridge 


5 


SC 


(«) 


Robbins 


Stanley 


6 


TC 


< K r) 


Proctor 


Proctor 


7 


SSC 


(k,t) 


Proctor 


Proctor 


8 


CSTC 


(kt, p) 


Robbins 


Mills- Robbins- Rumsey 


9 


esse 


(k,p) 


Robbins 


Kuperberg 


10 


TSSC 


(k, r, p) 


Robbins 


Andrews 



Table 2: Conjecturers and provers of plane partition formulas 
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Stanley 


Stembridge 


Group 


First q-formulator 


First q-enumerator 


1 


P 


(e) 


MacMahon 


MacMahon 


2 


S 


(r) 


MacMahon 


Andrews, Macdonald 


2' 


S 


(r) 


Bender-Knuth, Gordon 


Andrews , Gordon 


3 


cs 


(P) 


Macdonald 


Mills- Robbins-Rumsey 


4' 


TS 


(P,r) 


Andrews, Robbins 


still open 



Table 3: g-Conjecturers and g-provers of plane partition formulas 
a picture" of the partition, together with the back three sides of the box: 



We decree that the hexagon has angles of 27r/3 and congruent opposite sides of lengths a, b, and c (we will 
call such a hexagon semi- regular), and the rhombuses have angles of 7r/3 and 27r/3 and unit edge length 
(we will call such a rhombus a lozenge). Conversely, we can "look at" a tiling and "see" the corresponding 
partition. We invite the reader to arrive at this visual proof by staring at some representative rhombus 
tilings. 

The three operations r, p, and k have a simple interpretation in terms of tilings. The operation r is a 
reflection about a diagonal between two opposite vertices, p is a rotation by 27r/3 (clockwise by convention), 
and k is a rotation by it. For example, the number iV 9 (2a, 2a, 2a) is also the number of lozenge tilings of a 
regular hexagon which are invariant under rotation by 7r/3. 

Suppose we are given a tiling of a semi-regular hexagon H by lozenges. There is a unique tiling of H 
by unit equilateral triangles. We may form a graph Z whose vertices are these triangles, with two vertices 
connected when the triangles are adjacent. Each tiling by rhombuses corresponds to a (perfect) matching in 
Z. This brings us to the second idea, the permanent-determinant method, which can be used to count the 
number of matchings in a bipartite planar graph, and the Hafnian-Pfaffian method, a generalization for non- 
bipartite graphs. The general method is due to Kasteleyn [^0) B, but was motivated by an important special 
case independently discovered by Kasteleyn J| and Temperley and Fisher [fL5| , which in turn expanded on 
ideas of Hurst and Green 0] , Caianiello Q , and Kac and Ward . The permanent-determinant case of the 
method was noted by Percus jll] . 

We can use the permanent-determinant method in cases 1, 3, 6, and 8. The Hafnian-Pfafhan method 
applies in the remaining cases. 

This paper is a result of a collaboration with James Propp, who has maintained a steady interest in its 
contents and its completion, and the ideas presented here grew out of those in [Q, |5j- In particular I would 
like to thank him for correcting many minor errors in the final version. I would also like to thank William 
Jockusch, Richard Stanley, John Stembridge, and my former advisor, Andrew Casson, for encouragement 
and helpful discussions. 
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2 Definition of the general methods 



We begin with a brief overview of four functions of matrices known as the permanent, determinant, Hafnian, 
and Pfafhan. 

Let A be an n x n matrix. The permanent of A is a sum over permutations of n letters: 

n 

Per(A)= J2 Il^W 

aeS„ i=l 

The determinant is a signed sum of the same type: 

n 

Det(A)= ^(-irn^-W 

crGS„ 8=1 

In general, the permanent of a matrix is hard to compute (see Jl7|) and satisfies relatively few interest- 
ing identities, while the determinant is a fundamental object of all mathematics and can be computed in 
polynomial time. 

A (perfect) matching of a set S of 2n objects is a set of n disjoint pairs of elements of S. A matching is 
ordered or unordered depending on whether its members are ordered or unordered pairs. Let U n be the set 
of unordered matchings among 2n objects, and let O n be the set of ordered matchings. The Hafnian and 
the Pfaffian are polynomials which bear the same relationship to unordered matchings that the permanent 
and determinant do to ordered matchings, although the Pfaffian is most conveniently defined by ordered 
matchings. (Both the term "Hafnian" and the polynomial it denotes were devised by Caianiello, who explains 
that the Hafnian, like the element Hafnium, is named after the city of Copenhagen.) Like a permutation, 
an ordered matching has a sign: The matching (1, 2), (3, 4), ... , (2n — 1, 2n) is defined to be positive, and if 
er is a permutation of 2n elements and m is an ordered matching, then the sign of a(m) is defined to be the 
product of the signs of m and a. 

If A is a symmetric 2n x 2n matrix, the Hafnian of A is defined to be the sum: 

Hf(A)= yi n a ^ 

m£U n {ij}£ra 

If A is an antisymmetric 2n x 2n matrix, the Pfaffian of A is defined to be the sum: 

pf (^) 4^ { - i)m n a ^ 

After close scrutiny of this formula, it becomes clear that the use of oriented matchings is a notational device 
and the Pfaffian is really a polynomial with one term (with a leading coefficient of ±1) for each unoriented 
matching. The Hafnian is a generalization of the permanent and the Pfaffian is a generalization of the 
determinant, because if 



and 

A' = 






B 


B 








B 


—B 






then 

Hf {A) = Per(B) 

and 

Pf(A') = Det(B) 

Being a generalization of the permanent, the Hafnian is also in general intractible. But the Pfaffian satisfies 
the important identities 

Pf(5 T AB) = Det(B)Pf(A) 
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and 

Pf(A) 2 =Det(A), 

which are closely related to each other. For a proof of these assertions about the Pfaffian, see [H. It is also an 
instructive exercise to develop both the above definition of the Pfaffian and its properties from the following 
alternative definition: If the anti-symmetric matrix A is viewed as the matrix of a 2-form u>a e P 2 ", and 
\i is the standard volume form on M 2 ™, then: 

u>T = Pf(A)/x. 

Let Z be a bipartite graph with black vertices b±,...,b n and white vertices w\ . . . , w n . We define the 
bipartite adjacency matrix Mz of Z by setting (Mz)i,j to be the number of edges from bi to Wj. We see 
that the permanent of Mz is simply the number of matchings in Z. More generally, if Z w is a weighted 
version of Z, we define (Mz w )i,j to be the total weight of all edges from bi to Wj. The permanent of Mz w is 
the total weight of all matchings in Z w , where the weight of a matching in Z w is the product of the weights 
of the edges of the matching. 

Suppose that Z is planar and is given with a specific planar embedding, and suppose Z w is a weighted 
version of Z. We wish to obtain a new wighting Z w i, by changing the sign of some of the weights, so that: 

Det(M^,) =±Per(M z J (1) 

on a term-by-term basis. We define the weights of Z w to be the product of those of Z w and the weights of 
a signed graph Z±. We say that Z± is flat if for every face of Z (in the given embedding) that has 4fc sides, 
an odd number of sides are negative in Z±\ and for every face that has 4fc + 2 sides, an even number of sides 
negative in Z±. 

For the Hafnian-Pfaffian case, let Z be a weighted graph which is not necessarily bipartite. If Sz is the 
usual symmetric, weighted adjacency matrix of Z, then Hf(Sz) is the weighted sum of all matchings in Z. 

To have a Pfaffian, we need an anti-symmetric matrix. We let Z w be an oriented, weighted graph with 
vertices v±, . . . , v n . We define the anti-symmetric incidence matrix Az w of Z w by Aij = Wij — Wji, where 
Wij is the total weight of all edges from m to vj. 

Suppose that Z is an unoriented weighted graph with a given planar embedding. We wish to orient the 
edges of Z to obtain a graph Z w so that: 

Pf(A z -J = ±m(S z J (2) 

Once again, we use a consistency rule: We say that the orientation of Z w is flat if each face F of Z w has an 
odd number of edges which point in the clockwise direction around F. 

Kasteleyn's theorem states that every planar graph has a flat orientation, and that equation ^ holds if Z w 
is flat. These results have an excellent treatment in fllO[ which would not be improved if it were repeated here. 
However, just as the determinant is an interesting special case of the Pfaffian, the permanent-determinant 
method is an interesting special case of the Hafnian-Pfaffian method that admits a self-contained treatment 
that we present here: 

Theorem 2.1. Every planar, bipartite graph Z with an even number of vertices has at least one flat weighting 

z±. 

Proof. We view Z as embedded in the sphere and we will achieve flatness on the outside face as well as the 
others. (In fact, it is impossible to have flatness on all but one face on the sphere.) Let /o be the number 
of faces with 4fc sides, let fi be the number of faces with 4fc + 2 sides, and let e and v be the number of 
edges and vertices. Then e = fi mod 2 while v = mod 2. Therefore the Euler characteristic equation 
/o + /2 — e + v — 2 reduces to f% = mod 2, or in other words, an even number of faces have 4fc + 2 
sides. Choose a matching among these faces and connect their centers with paths on the sphere which are 
transverse to Z. If an edge of Z is crossed by n paths, give it a weight of (— l) n in Z±. Since every path 
enters and leaves a face with 4k sides the same number of times, Z± is flat at such a face. For a face with 
Ak + 2 sides, the number of entries has the opposite parity from the number of exits, therefore Z± is flat 
here also. □ 
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Theorem 2.2. If Z w is an arbitrary weighting of a planar bipartite graph Z , Z± is a flat signing, and Z w > 
is the product weighting, then: 

Det(Af Ziu/ ) = ±Per(M z J 

Proof. We would like to prove that all terms in Det(Af^ ,) have the same relative sign compared to the 
corresponding terms in Per(Mz TO ). It suffices to show that all terms in Det(Mz ± ) have the same sign. For 
simplicity we assume that Z has at most one edge connecting any pair of vertices. 

We first note that if a cycle c in Z, not necessarily the boundary of a face, encircles an even number of 
vertices, then the cycle has an even number negative edges if it has 4fc + 2 sides and an odd number if it 
has 4fc sides. (Proof: We define Z c> ± by removing all edges and vertices of Z± outside of c. We know that 
Z c .± is flat at every internal face and we want to show it is flat on the outside face as well. To do this we 
show that any sign ed g raph G± with an even number of vertices has an even number of non-flat sides. By 



the proof of lemma 2.1, we know that G± has an even number of faces with 4fc sides. If G± has the trivial 
signing, these are precisely the sides which are not flat. But if the sign of an edge is reversed, the flatness of 
precisely two sides is reversed.) 

Let mi and 7712 be two matchings in Z, and let o\ and 02 be the corresponding permutations. The set 
tyl\ U 7712 — mi H 7712 is a union of disjoint cycles, each encircling an even number of vertices of Z, and there 
is a sequence of matchings in Z connecting mi to 7712 such that two consecutive matchings in the sequence 
differ by only one cycle. We want to show that the terms for all of these matchings have the same sign. 
Therefore we can assume by induction that mi U 7712 — ttii PI 7712 consists of a single cycle c. 

Suppose that c has length 2n. (By construction it alternates between edges of 7711 and 7712 and therefore 
has even length.) Then the permutation a^ x a\ is an 71-cycle, which is an even permutation if 71 is odd and 
vice-versa. Therefore the relative signs of the permutations is exactly canceled by the relative signs of the 
edges of mi and 7712, which is given by the number of negative edges in c. □ 

The rows and columns of the matrices in this paper will in general be indexed by an unordered index set 
rather than by the integers from 1 to n for some 71. The Pfaffian of such a matrix is only defined up to sign, 
and since the rows and the columns may have different index sets, the determinant has the same ambiguity. 
Therefore we will work instead with the absolute determinant or absolute Pfaffian (the absolute value of the 
determinant of Pfaffian), which are still unambiguous quantities in this case. 



3 Enumeration in the simplest case 

In this section we give a complete proof of the formula for the number of plane partitions in case 1. The 
analysis of case 1 will set a pattern for the other cases and will also be used as a lemma. 

We will explicitly define matrices whose determinants or Pfaffians are the answers to the ten counting 
problems. 

Let H(a, b, c) denote a semi-regular hexagon with edge lengths a, 6, and c. Let T(a, b, c) be the tiling of 
H (o, b, c) by equilateral triangles of unit edge length. We can divide the triangles of T(a, b, c) into two kinds, 
V triangles and A triangles: 




G 



We define a graph Z(a, b, c) with vertex set T(a, b, c), where two triangles are connected by an edge if and 
only if they are adjacent. Z(a, b, c) might look like this: 



We give Z(a,b,c) the "default" weighting, whereby the weight of each edge is 1. Evidently, Z(a,b,c) is 
bipartite, and we can take the labels V and A as a standard bicoloring. 

We see that a matching in Z(a,b,c) corresponds to a plane partition in B(a,b,c). Thus, N\(a,b,c) is 
simply the number of matchings in Z(a, b, c). Let M(a, b, c) be the bipartite adjacency matrix of Z(a, b, c). 
We decree that the V triangles are the rows of M(a, b, c) and the A triangles are the columns. To apply the 
permanent-determinant method, we need to find a new matrix whose determinant equals Per(M(a, b, c)). 
In this case, all of the faces of Z(a, b, c) except the outside face have six sides. Therefore not changing the 
weights of Z(a, b, c) at all is a flat sign rule. In other words, 

|Det(M(o,6,c))| = Per(M(a, b, c)) = 7Vi(a,6,c). (3) 

We have obtained a determinant which we must evaluate. We proceed by row- and column-reduction 
and by induction on the area of H(a 1 b, c). 

We embed H(a, b, c) in H(a + 1, b + 1, c — 1) as follows: 




We embed T(a, b, c) in T(a + 1, b+ 1, c— 1) and M(a, b, c) in M(a + 1, b+ 1, c— 1) in the corresponding way. 
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We label the triangles in and adjacent to T(a + 1, b + 1, c — 1) — T(a, b, c) with the labels Xi, j/i, z,, and i: 




We can use each label twice, once for a V triangle and once for a A triangle. We use these labels to write 
out some pieces of the matrix M(a + 1, b + 1, c — 1). For example, let D be the submatrix whose rows are 
I and the x,'s and whose columns are the t/j's and Zj's. Note that D T also appears as a submatrix in the 
symmetrical position. In the case c = 4, the matrix D looks like this: 





2/1 


2/2 


2/3 


Z<3 


Zl 


Z2 


Z-3 


I 


1 








1 











Xl 


1 


1 








1 













1 


1 








1 





2^3 








1 











1 



If we multiply D by the matrix E: 





I 


Xl 


X2 


2^3 


I 


1 


-1 


1 


-1 


Xl 





1 


-1 


1 


X2 








1 


-1 


x 3 











1 



on the left, we obtain the product: 





yi 


2/2 


2/3 


zo 


Zl 


Z2 


Z3 


I 











1 


-1 


1 


-1 


Xl 


1 











1 


-1 


1 


X2 





1 











1 


-1 


X3 








1 











1 



Since E has determinant 1, and since the x, and I rows and columns are zero outside of D and D T , we can 
replace D and D T by ED and D T E T in M(a + 1, & + 1, c — 1) to obtain a new matrix M'(a, b, c) with the 
same determinant. 

The Xi and yi rows and columns are irrelevant to the absolute determinant of the matrix M'(a,b,c), 
because the Xi row (column) is the only row (column) in which the yi column (row) is non-zero, and the 
entries at the intersections of these rows and columns are all ±1. If we delete these rows and columns we 
obtain a matrix R(a, 6, c) which still has the same absolute determinant. This matrix is simply M{a, b, c) with 
an extra row and column labelled I. To complete the description of R(a,b,c), we note that i?(a, 6, c); jZi = 
R(a, b, c) Zi j = (— l) 1 and R(a, b, c)i_ t = or R(a, b, c) t j = where t is any triangle which is not a z,- L triangle. 

We are left with the problem of reducing the I row. Since M(a, b, c) is invertible (because, by induction 
its determinant is the number of plane partitions and there is at least one of them), we know that there is 
some linear combination of the M(a,b,c) rows which matches the I row exactly except for the R(a,b,c)i j i 
entry. We choose a coefficient a t for each V triangle t 6 T(a, b, c) so that: 

^2 a t R{a, b, c) t , u + R{a, b, c)j, u = 
t 
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for each A triangle u e T(a, b, c). After adding a t times the t row to the I row for each t, we obtain a new 
matrix R'(a, b, c) for which: 

R'(a, b, c\i = ^2 a tR( a , b, c) M 
t 

The matrix R'(a, b, c) still has the same absolute determinant as M(o + 1, b + 1, c — 1) and has the form: 





• \ 


M(a,b, c) 




... 


r ) 



where r — R'(a, b, c)ij. Therefore: 

|Det(M(a+l,&+l,c-l))| = \&(a,b,c)i,t\\Det(M(a,b,c))\, (4) 

which will prove the inductive step. For convenience, we define on — a Zi and fa — R(a, b, c)i tZi — R(a, b, c) Zit i, 
where i is an integer. Observe that fa = (— l) 1 . 

Let s,t,u, and v be four adjacent triangles in the middle of H(a, b, c) arranged as follows: 



s A t 

, / v\ 



We see that the corresponding a coefficients satisfy a symmetric Pascal triangle rule: 

a s + a t + a u = -R(a, b, c);,„ = 
The fa coefficients may be incorporated into this rule, because if s,t, and z% are arranged as follows: 



then: 

a s + a t + R{a, b, c)i lZi = a s + a t + fa = 

If we put each coefficient a t inside the triangle t, and each Pi under the A triangle Zi, we see that each 
number is the negative of the sum of the two above it: 




4 A. -4 A 4 / 
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(In this figure the coefficients have all been multiplied by 4 to clear denominators.) Thus, we get a set of 
linear equations for the a,'s in terms of /Vs: 

Strictly speaking, the summation over j should be from to c — 1 and the equation is only valid for i from 
to c — 1. However, if we declare that a.; = and (3i is an indeterminate for i < and i > c, then the 
equation will be valid for all z and the summation can be extended to go from — oo to oo. (Recall that 
out-of-bounds binomial coefficients are zero by common convention if the top index is a non-negative integer 
and the bottom index is an integer. Among other things this implies that the /3,'s must have finite support 
if the Qfi's do.) All subsequent summations will run over all of Z unless otherwise noted. 
We also get an expression for the final answer in terms of the ai and Pi coefficients: 



R'{a,b, c)i t i = 22, a iPi- 



The author originally obtained the solution to these equations with the aid of MACSYMA™, but with 
hindsight we may use generating functions to make the computation look easy. We will work in the formal 
power series ring M.[[x, y]][l/x, l/y] (which should not be confused with R((x, y)))- Keeping in mind that the 
ai and j3i are implicitly functions of a, b, and c, we define these generating functions: 



A(a,b, c;x) = a.jX l 

i 

oo 

A(a,b;x,y) = A(a, b, c; x)y°' 1 



c=0 



M{a,b-x) = {-l) a+b x- b (l + x) c 
M(a, b; x, y) = M(a, b; x) = (-1)°+^- V"" + vT +b 
B(a,b,c;x) = ^ (3 c -x-jX % 

i 

oo 

B(a,b:x,y) = ^B{a,b,c;x)y c ^ 1 



B(a, b, c; x) — M{a, b\ x)A(a, b, c; x), 



Then the equation 
which is equivalent to 

B(a, b; x, y) = M(a, b; x, y)A(a, b; x, y), 
is also equivalent to equation |^. By the definition of the j3i coefficients, 

1 



B(a,b;x,y) 



(l-xy)(l + y) 



where the ellipses here and below represent indeterminate or irrelevant terms with cither more factors of x 
than y or negative powers of x, and A(a, 6; x, y) is to be determined. Proceeding directly to the answer, we 
set: 

(_!)«+& 



A(a,b;x,y) 
We apply the binomial expansion theorem: 
M(a,b;x,y)A(a,b;x,y) 



(l-xyr+^l + y)b+i( a+ a b ) 

((l+y)-(l-xy)r+ b 
x b y a + b (l-xy) a + 1 (l + y) b + 1 ( a + b ) 

a+b ffl(-l)fr(l-, Ty )fc( 1+y) a+^ 

^ a;V +b (i - xy) a+1 (i + y) b+1 ( a V) 
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A given term in this sum with k < a, when expanded in x and y, has more factors of x than y, while a term 
with k > a has only negative powers of x. Therefore, only the k = a term contributes to the coefficients of 
B(a,b;x,y) that we have constrained. We compute: 

(-l) a 

M(a,b;x,y)A(a,b;x,y) = .. 



x b y a+b (l - xy)(l + y) 
1 



(l-xy)(l+y) 
= B(a,b;x,y), 

as desired. In the second line of the equation we again discard terms with more factors of x than y or 
negative powers of x. 
Finally, we recall: 

R'(a, b, c) M = cafti = A(a, b, c; -1) 

i 

and we compute: 

°° (_l)a+b 

V A(a, b, c; -l)y c_1 = A(a, b: -1, y) = '- 

Using the power series expansion: 



(1 - x) k 
we obtain: 



or equivalently, by equations || and |[ 

i\Ti(o + 1,6+ l,c- 1) /a + 6 + c\ If a + b 



Nx{a, 6,c) V c^ 1 

which is the desired ratio. 

For comparison with the other cases of enumeration, we introduce an extra subscript for many of the 
quantities defined here. We will use a\ i for oti, (3\^ for /3j, A\ for A, and i?i for _B. 



4 Determinants and Pfamans for matchings with symmetries 

In this section, we will analyze the set of matchings of a graph which are invariant under some group action. 
We will then apply this analysis to the particular group actions on the graph Z(a,b, c) under consideration. 
The result for each group will be a planar graph whose matchings correspond to plane partitions with the 
equivalent symmetry. We can then use the permanent-determinant or Hafnian-Pfaffian method to reduce 
each of the ten enumerations to questions about matrices. 

Suppose that Z is a finite graph with a distinguished vertex called "bachelorhood" . We define a matching 
with bachelors in Z to be a matching in which every vertex is matched to exactly one of its neighbors, except 
for the bachelorhood vertex, which can be matched to any subset of its neighbors. (Note that there may be 
more than one edge between the bachelorhood vertex and another vertex. We define the bachelor's degree 
of a vertex to be the number of edges connecting it to the bachelorhood vertex.) The concept of matchings 
with bachelors is a generalization of ordinary matchings in the sense that if Z is an ordinary graph, we let 
Z' be Z together with an isolated bachelorhood vertex, and a matching in Z will correspond to a matching 
with bachelors in Z' . 

Given a graph Z with a bachelorhood vertex, we wish to construct a graph Z' such that ordinary match- 
ings in Z' correspond to matchings with bachelors in Z. For this purpose we can replace the bachelorhood 
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vertex of Z by a subgraph which consists of either a row of triangles: 




bachelorhood 



or a row of triangles and an edge: 




bachelorhood 



The first subgraph has the property that if an odd number of vertices are matched to the outside, the 
remainder can be matched among themselves in a unique way, and the second subgraph has the analogous 
property for an even number of vertices. (The proof of this fact is left as an exercise to the reader.) Thus, we 
replace the bachelorhood vertex by one of these two graphs in such a way that the total number of vertices 
is even and call the result Z' . The graph Z' has exactly one matching for each matching with bachelors of 
Z. Moreover, if Z is a planar graph, it possible to construct Z' to be planar also. 

Suppose that Z is a graph with a bachelorhood vertex (which may be isolated) and a group G acts on Z . 
That is, each element of G is some permutation of the vertices that fixes the bachelorhood vertex and edges 
of Z that preserves the relation of an edge containing a vertex. We wish to construct a modified quotient 
graph Zj /G whose matchings are bijective with G-invariant matchings of Z . If e is an edge reversed by G, 
i.e. there exists an element of G which fixes e but switches its endpoints, we replace e by an edge from each 
of its endpoints to the bachelorhood vertex. We perform this operation for all edges reversed by G and call 
the result Z' . The action of G on Z extends to an action on Z', and the G-invariant matchings of Z' can 
be identified with the G-invariant matchings of Z. Suppose that e is an edge of Z' with a vertex v which is 
not the bachelorhood vertex such that the stabilizer of e does not contain the stabilizer of v. Then e cannot 
be a member of any G-invariant matching of Z' . Let Z" be Z' with all such edges removed. Let Z" ' /G be 
the usual quotient of a graph by a group action: The vertices of Z"/G are the G-orbits of vertices in Z" , 
and two vertices are connected by an edge if the corresponding orbits of vertices are connected by an orbit 
of edges. Finally, we define Z//G to be Z" /G. The graph Z//G has the property that: 

Lemma 4.1. The matchings of Z//G are in natural bisection with the G-invariant matchings of Z . 

The proof of this lemma is sketched by the definition above. 

It is not true in general that if Z is planar then Z//G is also. Indeed, if Z is the edge graph of a hexagonal 
prism: 




then Z has a fixed-point free involution such that the quotient is -^3,3, namely the involution consisting 
of rotation by 180 degrees and transposition of the inner and outer hexagons. However, suppose that 
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Z = Z(a,b,c) for some a, b, and c, with an isolated bachelorhood vertex added, and suppose that G = Gi 
for some 2 < i < 10. Then Z(a,b,c)/ /Gi is indeed a planar graph in these cases. We further modify Z/ '/G 
to remove the bachelorhood vertex and call the result Zi(a, b, c). If i is 3 or 9, or if i is 5 and a, b, and c all 
have the same parity, then Gi acts freely on Z(a, b, c) and the quotient graph is particularly simple. When 
kt G Gi, namely in cases 6,7,8, and 10, the vertices and edges of Z{a, b, c) fixed by kt or its conjugates map 
to a collection of disjoint edges in Zi(a,b,c), which can be ignored since there is only one matching in a 
graph which is itself a matching. When r G Gi, namely in cases 2, 4, 7, and 10, Zi(a,b,c) has a non-trivial 
bachelorhood vertex, so it is necessary to modify it in the manner described above to obtain a planar graph 
without a bachelorhood vertex. Finally, in cases 3, 6, and 8, the end result is a bipartite graph, and therefore 
the permanent-determinant method is applicable. In the other six cases, the Hafnian-Pfaffian method must 
be used. 



5 Determinants and Pfaffians for q- enumeration 

The previous two sections describe how to construct matrices whose determinants and Pfaffians enumerate 
plane partitions in each of the ten symmetry classes. The permanent-determinant method applied to weighted 
graphs yields matrices whose determinants and Pfaffians q-enumerate plane partitions in various ways. In 
particular, the natural (/-enumerations described in fll3| for cases 1,2,3, and 4 can all be expressed this way. 
In this section we give a determinant for (/-enumeration in the simplest case, case 1, as an illustration of the 
general principle. 

We define the g-weight of a plane partition to be q n if the plane partition has n elements. We wish to 
find the total q-weight of all plane partitions in an a x b x c box. We say that two plane partitions in such 
a box differ by an elementary move if they are the same except for one cube. Observe that we can connect 
any two plane partitions in a box by a sequence of elementary moves. If two plane partitions differ by an 
elementary move, the corresponding matchings in Z(a, b, c) differ on only one hexagon: 




In fact, there are two kinds of elementary moves: A raising elementary move is in which a cube is added, 
while a lowering elementary move is one in which a cube is deleted. 

We wish to weight the edges of Z{a, b, c) in such a way that the weight of a matching agrees with the 
weight of the corresponding plane partition. If we find such a weighting with the property that if we perform 
a raising elementary move on a matching, the weight goes up by a factor of q, then the total weight of all 
matchings will equal the total weight of all plane partitions up to a constant factor. That constant factor is 
the weight of the matching corresponding to the empty plane partition, because the q-weight of the empty 
plane partition is 1. To accomplish this we give all slanted edges in Z(a, b, c) and all of the left-most vertical 
edges a weight of 1. The weight of any other vertical edge is q times the weight of the edge immediately to 
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the left of it. The result is of the following form: 




The bipartite adjacency matrix of this graph has the desired determinant up to a normalization factor. 

The total q- weight is also an interesting statistic in cases 2, 3, and 4, and in each of these cases the same 
method yields a matrix whose determinant or Pfafnan is the desired statistic. However, there is no known 
formula for the total q- weight of all plane partitions in case 4 in this sense. For the purpose of understanding 
case 4, there is an alternative weighting scheme: The g-weight of a plane partition is defined to be q n if the 
plane partition contains n orbits under the required symmetry group. The total orbit-g-weight satisfies a 
nice formula in cases 2 and 4, and again in those two cases as well as in case 3 the above method yields 
a determinant or Pfaffian. The orbit-g-enumeration of plane partitions in case 4 remains the only open 
enumeration in the standard conjectures about symmetric plane partitions. 

6 Enumeration in a previously open case 

The ultimate goal of this section is to prove the formulas for N 3 (a,a,a), N 5 (2a,2b,2c), and Ng(2a, 2a, 2a) 
by counting matchings in Z$(a,a, a), Z§(2a,2b,2c), and Zg(2a,2a,2a). The hard part of the proof will be 
row- and column-reduction. The proof in case 1 will be used as a template and as a lemma. In addition, the 
three cases treated here together with case 1 form an analogy: 

case 1 is to case 3 
is to as is to 
case 5 is to case 9 

Recall that the vertices of Z(a, b, c) are the triangles in a tiling T(a, b, c) of a hexagon H(a, b, c): 




Since Z 3 (a,a,a) is a quotient graph of Z(a,a,a), its vertices are p-orbits {t, p(t), p 2 (t)} of elements of 
T(a, a, a). (One of the p-orbits in the picture is shaded.) Moreover, two vertices of Z${a, a, a) are connected 
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by an edge if and only if they are adjacent as p-orbits. (Note that there is a pair of p-orbits in the middle 
which are connected by two edges.) Observe that elements of a p-orbit are either all Vtriangles or Atriangles, 
so we can speak of V p-orbits and A p-orbits. In other words, the bicoloring of Z(a, a, a) carries over to a 
bicoloring of Z 3 (a, a, a). As in case 1, the trivial sign rule is flat, so if we let M 3 (a,a,a) be the bipartite 
adjacency matrix of Z 3 (a,a,a), its determinant is N 3 (a,a,a). We decree that the V p-orbits should index 
the rows of M 3 (a, a, a) and the A p-orbits should index the columns. 

The vertices of Zc,(2a,2b, 2c) are pairs of triangles {t,n(t)}, where one triangle is a V triangle and the 
other is a Atriangle. Since the graph Z 5 (2a, 2b, 2c) is not bipartite, we will need to examine its antisymmetric 
adjacency matrix under some flat orientation rule, and such a matrix M 5 (2a, 2b, 2c) would have one row and 
one column for every vertex of Z^(2a,2b,2c). Instead of having the rows and columns indexed by these 
vertices, we can as in case 1 have the rows of M 5 (2a,2b,2c) indexed by V triangles and the columns by 
A triangles, because each V triangle corresponds to a unique pair and vice- versa, and the same is true of the 
A triangles. With this indexing, we can take advantage of the fact that M 5 (2a, 2b, 2c) (with any orientation 
for Z 5 (2a,2b,2c)) is simultaneously the bipartite adjacency matrix for some (non-fiat) signed version of 
Z(2a, 2b, 2c). The signs must be chosen so that the an edge from t to u is the negative of the edge from nt to 
ku, and it must also correspond to a flat orientation. One such weighting is characterized by the following 
map: 




In this weighting, the edge between two adjacent triangles which are below the center of H(2a, 2b, 2c) has 
weight —1, as does the edge between two triangles which are to the right of the center if one is below the 
center and one is above the center. All other edges have weight 1. We let M 5 (2a, 2b, 2c) be the bipartite 
adjacency matrix of this weighted graph. It is antisymmetric and its Pfaffian is N§(2a, 2b, 2c). 

The same analysis holds for Zg(2a, 2a, 2a), except that the antisymmetric adjacency matrix of an oriented 
version of this graph is simultaneously the bipartite adjacency matrix of a weighted version of Z 3 (2a, 2a, 2a). 
We give Z 3 (2a,2a,2a) the following weighting: 




and we let Mg(2a, 2a, 2a) be the bipartite adjacency matrix of the result. One way to describe this weighting 
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is to say that if t and u are adjacent triangles, then: 

M 9 (2a, 2a, 2o){ t , P (t),p 2 (t)} ,{«, P («), P (« 2 )} = M 5 (2<z, 2a, 2a) t , u M 5 (2a, 2a, 2a) p( ^ p(u) M 5 (2a, 2a, 2a) p 2 (t)p 2 (u) 

The only exception to this rule is in the very center: If t and u are the two p-orbits of triangles in the very 
center of H(2a, 2a, 2a), then Mg(2a, 2a, 2a)t, u = 0. 

We embed M 3 (a, a, a) in M 3 (a+1, a+1, a+1), M 5 (2a, 26, 2c) in M 5 (2a+2, 26+2, 2c-2), and M 9 (2a, 2a, 2a) 
in M 9 (2a + 2, 2a + 2, 2a + 2). The corresponding embeddings of semi-regular hexagons look like this: 




Case 9 



In each case, we label many rows and columns of the larger matrix for the purpose of row-reduction by 
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labeling the corresponding triangles or p-orbits: 




Case 9 



With these labels, the row-reduction algorithm has the same basic plan as in case 1: For i = 3,5, or 
9, we define a matrix Ri{a,b,c) by multiplying the rows and columns, the I row and column, and (if it 
exists), the k row and column by a certain matrix with determinant 1. In the result the Xj and yj rows 
and columns will be irrelevant, and the only non-zero entries of the I and k rows will be in the Zj, I, and k 
columns. However, the actual values of those entries will be somewhat different. In all three cases, we define 
Pij — Ri(a, b, c)i_ Zj and in cases 5 and 9 we also define 7^ = Ri{a, b, c) Zj _k- 

In the final row- reduction step we will find a linear combination of the rows of Mi(a, b, c) which matches 
the I row exactly except for the R(a, b, c); ; ; coefficient and, if it is defined, the R(a, b, c)i^ coefficient. We 
choose a coefficient a^ t for each row t of Mi(a, b, c) so that: 

^ a lyt R(a, b, c)t, u + R(a, b, c)j,„ = 
t 

for each column u. In case 3, we add 0:3 j times the t row to the I for each t to obtain a matrix R' 3 (a, b, c). 
Just as in case 1, the result is: 

|Det(M 3 (o + 1, a + 1, a + 1))| = \R' 3 (a, a, a) M ||Det(M 3 (2a, 2a, 2a))| 

In cases 5 and 9, we add times the t row to the I row and then add a^t times the column to the I 
column to obtain R' 5 (2a, 2b, 2c) or R' g (2a, 2a, 2a). All operations performed to obtain these two matrices will 
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be symmetric with respect to rows and columns, and as a result these two matrices must be antisymmetric 
as M 5 (2a, 2b, 2c) and Mg(2a, 2a, 2a) are. The matrix R' 5 (2a, 2b, 2c) will have this form: 





• 




• 


M 5 (2a,2b,2c) 


• 




• 


... 


r 




-r 



where r = R' 5 (2a, 2b, 2c)i t k- From the form of the matrix we obtain a relation for the absolute Pfafhan: 

|Pf (M 5 (2o + 2, 2b + 2, 2c - 2))| = \R' 5 (2a, 2b, 2c)i, fc ||Pf (M 5 (2a, 2b, 2c))\ 

and by the same argument for case 9, 

|Pf(M 9 (2a + 2,2a + 2,2a + 2))| = \R' g (2a, 2a, 2a); !fc | |Pf (M 9 (2a, 2a, 2a))| 

To derive the actual entries of the matrices Ri(a, b, c) and i?-(a, b, c) and the coefficients ctjj = cti iZj , we 
proceed on a case-by-case basis. As before, once the a, ; j's are chosen, the values of the other a, ;t 's can be 
determined in a simple manner. 



6.1 Case 3 

In case 3, the submatrix indexed by Xi and / rows and yi and z\ columns looks the same as in case 1, but 
we must include the / and x a columns because M^(a + 1, a + 1, a + l)j,x arL d Ms(a + 1, a + 1, a + l)x„,i are 
non-zero. Here is the submatrix D which covers all non-zero entries of the I and Xi rows in the case a = 3: 





X3 


yi 


2/2 


2/3 


zo 


Z\ 


Z2 


I 


I 


1 


i 








1 











Xl 





i 


l 








1 








X2 








l 


1 








1 





%3 











1 











1 



If we multiply D by the matrix E: 





I 


Xl 


x 2 


2^3 


I 


1 


-1 


1 


-1 


Xl 





1 


-1 


1 


X 2 








1 


-1 


x 3 











1 



on the left, we obtain the product: 





2/i 


2/2 


2/3 


zo 


Zl 


Z2 


I 


I 











1 


-1 


1 


-1 


Xl 


1 











1 


-1 


1 


X 2 





1 











1 


-1 


X3 








1 











1 



Since E has determinant 1, and since the x^ and I rows are zero outside of D, we can replace D by 
ED in M^(a + 1, a + 1, a + 1) without altering the determinant. However, notice that in case 3 there is an 
intersection of a single entry, M^{a, a, a)i t i, between D and Dt, and multiplying D on the left by E alters this 
entry. Nevertheless we still multiply on the right by E T . We again delete the Xi and yi rows and columns 
to obtain i?3(a,a, a). This matrix consists of M3(a,a,a) plus the non-zero entries 03,i = i?3(a, a, a)z, Zi = 
J?3(a, a, a) Zi ,i = (— l) 4 and also Rs(a, a, a)i,i = (— l) a 2. We are left with the derivation of the Q3,i's. 
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We can interpret the a^j coefficients as a function on individual triangles in T(a, a, a) rather than p-orbits 
by the equation: 

a 3,t = a3,{t,p(t), P 2 (t)} 

Like the ai,t's, the a^jS satisfy the symmetric Pascal rule in H(a,a,a), but the boundary conditions are 
different. To extend the symmetric Pascal rule to the boundary, we must place the fo^ coefficients on all of 
the triangles in the p-orbits of Xi, not just on the bottom row: 




10 -10 10 



We obtain the solution by setting: 

&3,t = Ot\,t + a l,p(t) + a l,p 2 {t) 

for all triangles t G T(a, a, a). For integers i, we also define a[ i = a\ tP t Zi \ and a'{ i = a\ tP 2t Zi y 
We expand R' 3 (a,a,a)i t i: 

R' 3 (a, a, a)i,i = R 3 (a, a, a) M + ^(-1)^1,, + ^(-1) V M + 

i i i 

= (-ir2 + (-l) 3 - 1 ( a 3a i yQ+2^(-l)V M (7) 

by equation ^| and by the symmetry a.\, TK {t) = (— l) a+1 ai.t, which implies that a'{ i = {— l) a+1 a' 1 a _ 1 _ i - 
From the position of p(zi) in H(a, a, a) (recall that p is a clockwise rotation), we obtain the formula: 



i 



Therefore: 



2a - 1 - j ' J x ' ^ V 2a - j 

Recall that f5\-i is the x a coefficient of the polynomial Bi(a, a, a; x) and the x a y a ~ 1 coefficient of Bi(a, a; x, y). 
We compute: 

B 1 (a,a,a;x,y) = M 1 {a,a;x,y)A 1 [a,a,a;x,y) - 



x a y 2a (l - xy) a+1 (1 + y) a+1 ( 2 a °) 



v {^)(-\f(l-xf(\+y)^- k 

4^V a (i-^) a+1 (i + y) a+1 (a a ) 
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Both the k = term and the k — a term, but no other terms, contribute to the x a y a 1 coefficient: 

+ (-!)■ 



B 1 (a,a,a;x,y) 



a y 2a (l - xy) a + l ( 2 *) x a y 2a (l - xy)(l + y) 

a..a-l ( f 3a \ I ( 2a> 



= X ~ y ~ '>> a)/[a '- 1 



(8) 



Substituting equation ^| in equation [7], we have: 

3d \ / / 2a \ „ / 3a\ / /2a 



This is the desired ratio up to sign. 



3a + 2 / 3a \ / /2a 



2 _ 

a / /V a / / a \a — 1 / / \ a 



6.2 Case 5 

In the two PfafRan cases, the Xi triangles near a given edge of the hexagon form two rows rather than one 
and as a result the initial row-reduction step is more complicated. As before, we consider the submatrix D 
of Ms{2a + 2, 2b + 2, 2c — 2) consisting of the X{, I, and k rows and the yi and Zi columns. Here is D when 
c = 2: 





Vi 


2/2 


V3 


2/4 


2/5 


2/6 


2/7 


2/8 


2/9 




Zl 


Z2 


Z3 


Z4 


Z5 


k 


-1 











-1 
































Xi 


-1 


-1 











-1 





























%2 





-1 - 


-1 











-1 


























X3 





- 


-1 


-1 











-1 























X4 











-1 














-1 




















I 














-1 














-1 

















X 5 














-1 


-1 














-1 














x 6 

















-1 


-1 














-1 











X 7 




















-1 


-1 














-1 








X$ 























-1 


-1 














-1 





X 9 


























-1 

















-1 


Note that — D T also appears as a submatrix of A/5 (2a 


+ 2,2&+2, 


2c- 


2). If we multiply D 


of the form: 






































k 


Xl 


X2 


X3 


24 


I 


£5 


X6 


XI 


Xs 


Xg 










k 


1 


-1 


1 


-1 


1 





-1 


2 - 


-3 


4 


-5 










Xl 





1 


-1 


1 


-1 








-1 


2 


-3 


4 










Xi 








1 


-1 


1 








- 


-1 


2 


-3 










X3 











1 


-1 














-1 


2 










X4 














1 

















-1 










I 

















1 


-1 


1 - 


-1 


1 


-1 










x 5 




















1 


-1 


1 


-1 


1 










x 6 























1 - 


-1 


1 


-1 










X7 


























1 


-1 


1 










X 8 





























1 


-1 










Xg 
































1 







We obtain: 
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Vl 


J/2 


2/3 


J/4 


J/5 


ye 


y~i 


J/8 


J/9 


Z{) 


Zl 


Z2 


Z3 


Zi 


Zb 


k 
































1 


-2 


3 


-4 


5 


X\ 


-1 
































1 


-2 


3 


-4 


X2 





-1 
































1 


-2 


3 


xz 








-1 
































1 


-2 


X4 











-1 
































1 


I 





























-1 


1 


-1 


1 


-1 


1 


x 5 














-1 

















-1 


1 


-1 


1 


-1 


x 6 

















-1 

















-1 


1 


-1 


1 


X7 




















-1 

















-1 


1 


-1 


x$ 























-1 

















-1 


1 


Xg 


























-1 

















-1 



We replace D by ED and -D T by -D T E T in M 5 (2a + 2, 26 + 2, 2c- 2), and then delete the Xi and yi rows 
and columns to obtain R^(2a, 26, 2c), which still has the same absolute Pfaffian. In addition to M 5 (2a, 2b, 2c), 
the matrix R${2a, 2b, 2c) has the non-zero entries (5^^ = R^{2a, 2b, 2c)i^ Zi = —R^{2a, 2b, 2c) Zit i = — (— l) 1 and 
75 ,i = -R 5 (2a,2b,2c) k , Zi =i(-l) 1 . 

The 015,4 coefficients, as a function on the V triangles of T(2a, 2b, 2c), satisfy the symmetric Pascal rule 
everywhere except along a "branch cut" going from the center of T(2a, 2b, 2c) to the right edge, where they 
satisfy the ordinary Pascal rule: 




We label the row of triangles just below the branch cut as indicated and we define 015^ = 0-5^. 
We define the following generating functions: 

A 5 (2a, 26,2c; x) = ^Va 5jl 

i 

A 5 (2a, 26, 2c; x) = ^ x *«5,i 

i 

B 5 (2a,2b,2c;x) = ^V/^c-i 

i 

For convenience in dealing with the branch cut, we define the function 9(n) to be —1 when n > and 1 
when n < 0. We define the functions <r and a on Laurent polynomials by: 

a(P(x)) = Y.°^ xl 

i 

a(P(x)) = ]>>«^ 

i#0 
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where: 

/':<•:■ 

i 

We examine the equations that the a^ i coefficients satisfy: 



We can put equations |p and [l^ in generating function form: 

3 5 (2a,26, 2c; x) = cr(x~ c Afi(a, 6; x)A 5 (2a, 26, 2c; x)) 
5 5 (2a,26,2c;x) = x c Afi(a, 6; x)A 5 (2a, 26, 2c; x) 

As before, for i > 2c — 1 or i < 0, 015^ — 0, but fei is indeterminate. 

As usual, we jump to the solution (originally found with the aid of MATHEMATICA™) by setting 

A 5 (2a, 2b, 2c; x) = (1 + a:)Ai(a, 6, c; x) 2 . 

Before demonstrating that this is the solution, we define the polynomials Bl{o,, b, c; x) and B//(a, 6, c; x) by: 

Si (a, 6, c; x) = Mi (a, 6; x)Ai(a, 6, c; x) = (-l) c_1 — — h B L {a, b, c; x) + B H {a, 6, c; x), (11) 

1 + x 

where Bh(ci, 6, c; x) has no terms with exponents less than c, and B^(a, 6, c;x) consists only of terms with 
negative exponents. By analogy with B\{a, 6, c; x), we also define: 

B H (a,b,c;x) = }^ H ,c-i-jX l 

i 

B L (a,b,c;x) = ^ Pl,c-i^x 1 

i 

Using the fact that the terms of A\{a, b, c; x) have degrees from to c — 1 inclusive, we can expand 
A 5 (2a,26,2c): 

A 5 (2a, 26, 2c; x) 

= a(x~ c M 1 (a 1 6; x)(l + x)Ai(a, 6, c; x) 2 ) 

(-l) c -! 

= cr(x c (l + x)A x (a, 6, c;x)(— — ■ — - + S L (a, 6, c; x) + B#(a, 6, c; x))) 

(1 + x) 

(-l) c -! 

= x c (l +x)Ai(a, 6,c;x)( — — — - + B L (a,b,c;x) - B H (a,b,c;x)) 

(1+x) 

The result is a cancellation of the "error terms" Bh and £?£,: 

x c Mi(a, 6; x)A 5 (2a, 26, 2c; x) 

f-l) c_1 (-l) c_1 

= (1 + x)(— — ■ h B L (a,b,c;x) + B H (a,b,c;x))(—— h B L (a,b,c;x) - B H (a,b,c;x)) 

1+x 1+x 

= — j- — + (-l^^B^a, 6, c; x) + (1 + x)( J B L (a, 6, c; x) 2 - S ff (a, 6, c; x) 2 ) 
1 + x 

= B 5 (2a,2b,2c;x) (12) 



Finally, since 75^ = we see that: 

cL4 5 (2a,26, 2c; x) 



R' 5 (2a,2b, 2c) 



= A(a,b,q-l) 2 =R'(a,b,c)l l (13) 

x=—l 
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6.3 Case 9 



Case 9 is similar to case 3 in the sense that the positions in which the k, I, and xi rows are non-zero in the 
matrix M 9 (2a + 2, 2a + 2, 2a + 2) include the I and k columns as well as the j/j and z, columns. We let D be 
the submatrix consisting of all of these rows and columns. Here is D when a = 4: 





yi 


J/2 


J/3 


J/4 


J/5 


J/6 


J/7 


J/8 


J/9 


z 


Z\ 


Z2 


Z3 


I 


k 


k 


-l 











-1 
































X\ 


-l 


-1 











-1 





























X2 





-1 


-1 











-1 


























x 3 








-1 


-1 











-1 























X4 











1 














1 




















I 














-1 














-1 

















x 5 














-1 


-1 














-1 














x 6 

















-1 


-1 














-1 











X 7 




















-1 


-1 














-1 








X 8 























1 


1 














1 





Xg 


























1 

















1 



Some of the entries are positive, because unlike case 5, in case 9 the V triangles x ai X2 a , and X2 a +i he 
outside of the negative region in the sign rule map for the matrix M 9 (2a + 2, 2a + 2, 2a + 2). 

Note that — D T also appears as a submatrix of Mg(2a + 2,2a + 2,2a + 2) but that the two matrices 
intersect in the 2x2 matrix formed by I and k. We multiply D by a matrix E of the form: 





k 


Xl 


X2 


X3 


X4 


I 


x 5 


xe 


x 7 


x & 


Xg 


k 


1 


-1 


1 


-1 


-1 





-1 


2 


-3 


-4 


5 


Xl 





1 


-1 


1 


1 








-1 


2 


3 


-4 


X2 








1 


-1 


-1 











-1 


-2 


3 


X3 











1 


1 














1 


-2 


X4 














-1 



















I 

















1 


-1 


1 


-1 


-1 




X5 




















1 


-1 


1 


1 




Xq 























1 


-1 


-1 




X 7 


























1 


1 




x$ 





























-1 




xg 



































This matrix is the same as the one in case 5, except that the x a , x 2a , and X2a+i columns have been 
negated. The product in this case is: 





Vi 


yi 


2/3 


yi 


y5 


2/6 


y-i 


y& 


2/9 


zo 


Zl 


Z2 


Z3 


I 


k 


k 
































1 


-2 


3 


-4 


5 


Xl 


-1 
































1 


-2 


3 


-4 


X2 





-1 
































1 


-2 


3 


X3 








-1 
































1 


-2 


X4 











1 


































I 





























-1 


1 


-1 


1 


-1 




x 5 














-1 

















-1 


1 


-1 


1 




x 6 

















-1 

















-1 


1 


-1 




XI 




















-1 

















-1 


1 




x$ 























-1 

















-1 




Xg 


























-1 




















In the matrix Mg(2a + 2, 2a + 2, 2a + 2), we multiply D by E on the left, which affects the l-k submatrix 
of —D T , and then we multiply the modified version of — D T by E T on the right. After multiplying on the 
left, the l-k submatrix is in general: 
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k 


I 


k 


2a- 


h 1 


-2a 


I 




1 


-1 



and after multiplying on the right, it becomes: 







k 


I 


k 







-2a- 1 


I 


2a- 


hi 






We delete the Xi and iji rows and columns and call the result Rg(2a, 2a, 2a). In addition to the 
above submatrix and the submatrix Mg(2a,2a,2a), Rg(2a,2a,2a) also has the non-zero entries (3gj = 
R 9 {2a, 2b, 2c)i iZt = -Rg(2a, 2b, 2c) %a = -(-l) 4 and 79,, = -R 9 (2a, 2b, 2c) k , Zi = i(-l)\ 

As in case 3, we interpret the agj coefficients as a function on triangles in T(2a, 2a, 2a) by the following 
equation: 

ag,t = "9,(4, p(t), p 2(t)} 

The resulting coefficients satisfy the symmetric Pascal rule everywhere except along three branch cuts instead 
of one: 




We can recycle the computations of case 5 by setting 0:5^ = a^j for all t 6 T(2a, 2a, 2a) except those which 
lie in the region R between the two branch cuts on the left side. For t £ R, we set 0:5^ = — a§.t- We see 
that the edge's satisfy the same rule as the cist's, except with different boundary conditions. Following case 
3, we can obtain the former from the latter by symmetrization: 

otg,t = "5,t + "5,p(t) + «5,p(t) 

We define: 

a ki = «5,p(zi) = a 5,p(zi) 
a 5,i = ®5,p 2 (zi) = -a5, P 2 0,) 
(Recall that p is a clockwise rotation.) We expand R 9 (2a,2a,2a)i tk '- 

R' g (2a, 2a, 2a) Lk = R 9 (2a, 2a, 2a)i, k + ^ i(-l)V* + ^ i{-l) l a' 5 l + ^ i(-l) l < 4 

i i i 

= 2a + l+Q a ^y^ 2 +J2(-iy^-^-l)a' 5tl (14) 

from equation [l^ and from the symmetry a§.t — -05,^7(1)1 which holds for triangles t below the center of 
T(2a, 2a, 2a) and which implies that a'l i = — a' 5 2a —i-i- From the position of p(zi) in T(2a, 2a, 2a), we learn 
that: 

a 5,i = \L 



4a - 1 - j 



"5j 
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Therefore: 



^(-l) i (2i-2a-l)a' 5)i = £(2z - 2a - 1) ( 2 " + l _ ■ ) 



4a + 1 - j 

4a \ _/ 4a 
4a - j7 V 4a + 1 - 3 



£((2a-l) 



"5, 



= 2/3 5 ,_ 2 - (2a - l)/3 B ,-i (15) 

We want to know (j§ _i and P5.-2, which are the x 2a and cc 2a+1 coefficients of B§{2a, 2a, 2a; x). If we combine 
equations [ll] and [l^, we obtain: 

B 5 (2a, 2a, 2a; a;) = (1 + x)Bh (a, a, a; x) 2 + . . . 

1 + x 

= ...+x 2a - x 2a+1 + . . . - (1 + x)(j3 H ^ lX a + (3 H ,-2X a+1 + ...) 2 

= (l-[3 2 H ,_ 1 )x 2a + (-l-l3 H ,-i(2p H ^ 2 + f3 H ^ 1 ))x 2a+1 + ... (16) 

Putting together equations [l4[ [l5[ and [ll| we obtain: 

a 2 /3a 4 
(3a + l) 2 a 



3a+ 1\ 2 //2a> 2 



iZ&(2o,2o,2o) t ,fc= /Q „ , lV3 ( „ ) /(„) -j9H,-i(4^H,-a-(2a-3) J 9 Hi _i) 



We know from equation g that: 

^-eD/CD^r: 1 )/! 2 ; 

We can obtain @h,-2 hi the same way we derived Ph.-i- 

{1 + yf- 1 2a{l + y) a -2 (-1) Q 



Bx(a,a,a;x,y) 



Thus: 



and: 



xV a (l -xy) a+1 ( 2 Q a ) x a y 2a (l -a;y) a ( 2 a a ) "' zV a (l + y)(l - xy) 



n.2 



3a + 1 

i?g(a, a, a); !fe = 



(3a +1) 2 V a 
This is the desired answer. 





/C: 






2a^ 


i'-C 


a > 





9a 2 + 6a + 1 /3a + l\ 2 / /2a\ 2 /3a + l\ 2 / /2a^ 2 



a 
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